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This  report  describes  the  results  of  the  auroral  magnetic  aspect 
sensitivity  studies  being  performed  by  M&S  Computing,  Inc.  , under 
Contract  No.  DASG60-74-C-0026  for  the  Ballistic  Missile  Defense 
Systems  Command  in  Huntsville,  Alabama.  PAR  Aurc  ral  Study,  Volume 
V,  dated  October  19,  1976,  constitutes  M&S  Computing's  Report  No. 
76-0052. 


The  purpose  of  this  investigation  was  to  determine  the  relationship 
between  the  radar  reflectivity  of  the  aurora  borealis  and  the  angle  between 
the  magnetic  field  lines  and  the  radar  propagation  path.  The  mass  of  high 
resolution  data  collected  during  the  11/2  years  of  PAR  auroral  studies 
has  provided  the  means  to  study  this  relationship  and  achieve  results 
previously  unattainable. 
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1.  INTRODUCTION 


The  purpose  of  this  report  is  to  describe  the  measurements  of  the  rela- 
tionship of  the  apparent  radar  auroral  reflectivity  to  the  magnetic  aspect  angle, 
the  angle  between  the  propagation  path  and  the  magnetic  field  lines  in  the  vicinity 
of  the  apparent  source  of  the  reflection.  This  research  has  been  performed  as 
part  of  the  PAR  Auroral  Study  which  began  in  early  1975.  Data  collection  was 
performed  on  several  occasions  in  the  fall  of  1975  and  in  the  spring  of  1976. 

Over  200  reels  of  magnetic  tape  were  recorded  altogether,  representing  about 
40  hours  of  auroral  observation.  On  four  of  these  data  gathering  periods,  con- 
tinuous observation  and  recording  was  performed  throughout  the  night.  The 
recorded  data  consisted  of  auroral  backscatter  measurements  as  well  as  simul- 
taneous satellite  tracking  data.  In  addition  to  the  normal  data  collecting  activi- 
ties, a number  of  special  tests  were  performed  to  evaluate  the  effects  of  auroral 
clutter  on  ‘■he  sidelobe  blanker  and  to  measure  the  effects  of  active  auroral  noise 
generation  on  the  radar  noise  level. 

A number  of  special  computer  programs,  including  a versatile  graphics 
package,  were  developed  for  use  as  tools  in  this  analysis.  A more  detailed 
description  of  these  tools,  and  of  the  PAR  auroral  study  containing  numerous 
auroral  maps  and  special  plots,  has  already  been  published  in  three  volumes 
[References  1,  2,  and  3]  which  preceded  this  document.  The  development  of 
this  report  depends  heavily  on  background  information  supplied  in  those  volumes, 
which  are  identified  in  the  List  of  References.  Copies  of  these  documents  may 
be  obtained  from  the  Army  Ballistic  Missile  Defense  Systems  Command  in 
Huntsville,  Alabama. 

1.  1 Features  of  the  PAR 

The  PAR  is  particularly  well-suited  for  auroral  backscatter  studies. 
Among  its  features  are: 

• An  operating  frequency  in  the  450  MHz  band. 

• Phased  array  steering,  providing  rapid  wide-angle  scan 
capabilities. 

• Exceptionally  high-peak  and  average  transmit  powers,  allow- 
ing a range  coverage  well  over  2000  km. 

• Excellent  sensiti\.ty  and  noise  performance. 

• Narrow  receive  beam,  1.5°. 

• Excellent  range  resolution  of  better  than  1.5  km  accompanied 
by  very  high  transmit  pulse  energy  which  is  attained  by  pulse 
compression  techniques. 
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Optimum  face  orientation  to  minimize  distortion  and  steering 
loss  in  direction  of  maximum  auroral  backscatter.  Figure  1-1 
shows  the  location  of  the  PAR  along  with  its  geomagnetic  off- 
perpendicular  and  L-shell  contours  at  an  altitude  of  105  km. 

High  performance  data  processing  and  recording  capabilities. 


Capability  to  perform  simultaneous  measurements  of  auroral 
backscatter  and  satellite  tracking. 

Sidelobe  suppression  capabilities  to  minimize  unwanted  side- 
lobe  returns. 


Either  right  or  left  circular  polarization  which  can  be  selected 
for  both  transmit  and  receive. 


>lications  of  This  PAR  Stud\ 


With  these  capabilities,  auroral  data  from  the  PAR  provides  its  users 
with  a clear  picture  of  its  interaction  with  the  aurora  and  the  effect  of  the  aurora 
on  its  performance.  It  can  provide  valuable  engineering  data  to  the  designers 
of  future  radar  systems.  In  addition,  it  is  hoped  that  this  study  can  provide  the 
scientific  community  with  certain  valuable  geophysical  data  concerning  the 
aurora. 


1. 3 Overview 

Section  2 of  this  report  introduces  the  subject  of  aspect  sensitivity,  what 
it  is,  and  why  it  is  of  interest.  This  section  also  provides  a short  historical 
background  describing  the  results  of  some  earlier  aspect  studies. 

Section  3 describes  the  techniques  used  to  determine  the  magnetic  aspect 
relationships,  including  the  data  selection  process,  the  data  reduction  task,  and 
the  identification  and  handling  of  a number  of  analytical  problems. 


Section  4 presents  the  test  cases,  their  supporting  data,  and  the  resultant 
scatter  plots  which  show  the  aspect  response  curves. 

Section  5 describes  the  resuLts  of  the  analysis  and  summarizes  the  data 
from  the  test  cases.  Some  interesting  observations  have  been  made;  the  reader 
is  cautioned,  however,  that  some  of  these  observations  may  be  illusory.  This 
situation  is  discussed  in  further  detail  in  Section  6,  which  provides  possible 
explanations  for  the  observations. 
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ASPECT  STUDY  BACKGROUND 


The  phenomena  of  radio  wave  reflections  from  auroral  ionization  has 
been  observed  since  the  early  1930's.  These  reflections  have  been  measured 
throughout  the  UHF  range  of  30  MHz  to  3000  MHz,  but  most  of  the  measure- 
ments have  been  made  in  the  50-500  MHz  region.  At  frequencies  lower  than 
these,  the  radar  equipment  becomes  large,  cumbersome,  and  very  expensive. 

At  higher  frequencies,  the  reflectivity  of  the  auroral  irregularities  is  reduced, 
resulting  in  very  weak  echoes,  as  shown  in  Figure  2-1. 

One  particularly  significant  fact  about  auroral  backscatter  is  that  the 
intensity  of  the  observed  returns  depends  upon  the  angle  between  the  radar  beam 
and  the  magnetic  field  lines  at  the  reflecting  point.  This  phenomena,  called 
geomagnetic  aspect  dependence,  has  long  been  of  interest  to  scientists.  A num- 
ber of  theories  attempting  to  explain  the  geomagnetic  aspect  dependence  have 
been  developed  over  the  years.  Most  of  the  earlier  theories  have  now  been 
abandoned  in  favor  of  thjse  which  relate  the  irregular  plasma  structure  respon- 
sible for  the  observed  aspect  dependence  to  magnetohydrodynamic  wave  action. 

This  wave  action  is  thought  to  result  from  instabilities  in  the  plasma  which 
occupies  the  ionospheric  E region.  When  aurora  occurs,  the  plasma  density  is 
enhanced  and  also  subject  to  dynamic  interactions  related  to  the  geomagnetic 
disturbances  and  particle  precipitation.  The  resulting  wave  action  forms  the 
field-aligned  irregularities  which,  in  turn,  serve  as  reflectors  of  radar  signals. 

The  radar  signals  encounter  cyclic  changes  in  the  index  of  refraction  in  the  pro- 
pagation path  due  to  the  intensification  and  rarification  of  the  wave  fronts.  Basic 
electromagnetic  theory  predicts  that  energy  will  be  reflected  from  these  irregu- 
larities. The  echoes  detected  at  the  radar  consist  of  reflections  from  a great 
number  of  these  closely  spaced  irregularities. 

The  precise  processes  which  control  the  formation  and  dynamics  of  the 
magnetohydrodynamic  wave  action  are  not  completely  understood,  but  one  popular 
explanation  has  been  proposed,  known  as  "Two-Stream  Instability  " [5],  It  is 
not  within  the  scope  of  this  report  to  discuss  in  detail  the  various  theories  which 
attempt  to  explain  the  phenomena  but  primarily  to  provide  data  which  may  be 
used  to  support  and  advance  those  theories. 

Scientists  have,  for  a number  of  years,  been  trying  to  determine  the  geo- 
magnetic aspect  sensitivity  relationship  using  radar  measurements.  Most  of 
these  studies  were  hampered  by  the  equipment  used,  low  power,  broad  beam- 
width,  and  poor  range  resolution.  W.  G.  Chestnut,  a prominent  figure  in  auroral 
backscatter  studies,  made  aspect  measurements  at  numerous  frequencies.  An  ex- 
ample of  his  results  is  shown  in  Figure  2-2.  The  solid  line  represents  the  measured 
aspect  data,  and  the  dashed  line  represents  a prediction  based  on  the  measured 
data  of  what  the  aspect  relation  was  expected  to  be. 
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Preceding  page  blank 


A more  recent  effort  which  employed  considerably  more  advanced  radar 
equipment  was  undertaken  at  the  Prince  Albert  Radar  Laboratory  (PARL)  in 
Saskatchewan,  Canada.  This  study,  a joint  effort  of  Bell  Telephone  Laboratories, 
MIT  Lincoln  Laboratories,  and  Stanford  Research  Institute,  is  especially  signi- 
ficant because  of  the  similarities  between  the  PARL  and  the  PAR.  The  results 
of  this  study  have  been  published  in  two  documents,  [4]  and  [?].  One  significant 
factor  of  note  in  the  PARL  reports  is  the  use  of  "reflectivity"  to  describe  the 
volume  backscatter  properties  of  the  aurora.  The  use  of  reflectivity,  or  effec- 
tive radar  cross-section  per  unit  volume,  normalizes  the  range  and  pulse-width 
dependence  of  the  range  resolution  cell  as  well  as  most  of  the  constant  radar 
parameters.  This  technique  is  used  almost  exclusively  for  the  processing  of 
meteorological  radar  data  providing  those  researchers  with  an  improved  level 
of  standardization.  For  a discussion  of  the  techniques  employed  in  this  similar 
field,  see  [8]. 

Auroral  measurements,  described  in  terms  of  reflectivity,  allow  the 
direct  comparison  of  data  from  different  radars,  if  these  comparisons  are  made 
carefully  in  regard  to  the  magnetic  aspect  dependence.  Figure  2-3  is  an  example 
of  some  of  the  PARL  aspect  data  using  the  same  type  of  chirp  waveform  as  that 
used  by  the  PAR.  Another  set  of  aspect  response  curves  is  shown  in  Figure  2-4. 
These  curves  were  published  as  part  of  an  excellent  study  of  radar  auroral 
phenomena  [7], 

Unfortunately,  the  PARL  was  located  in  such  a way  that  measurements 
could  not  be  made  closer  than  about  4.5°  to  normal  magnetic  incidence.  Three 
significant  observations  were  made,  however.  The  first  was  that  the  "apparent" 
magnetic  aspect  sensitivity  depends  on  the  planetary  geomagnetic  disturbance 
index  (Kp).  This  index  indicates  changes  in  the  ionospheric  current  system. 

The  second  observation  was  that  the  slope  of  the  magnetic  aspect  sensitivity  tends 
to  become  progressively  flatter  (less  aspect- sensitive ) at  greater  off -perpendi- 
cular angles.  A third  observation  showed  that  the  reflectivity  depends  on  the 
effective  radar  pulse  width.  It  has  been  reasoned  that  this  last  phenomena  is 
merely  a manifestation  of  the  non-homogeneity  of  the  auroral  mass.  It  has  been 
suggested  [7]  that  with  the  longer  pulses,  the  range  resolution  cells  contain 
voids  and  weakly  reflecting  regions.  Voids  such  as  these  in  the  range  resolu- 
tion cell  decrease  the  average  reflectivity.  These  observations  are  all  supported 
by  the  PAR  Auroral  Study. 
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3.  DETERMINATION  OF  THE  MAGNETIC  ASPECT  RELATIONSHIP 


The  primary  intent  of  this  portion  of  the  study  was  to  determine  the  mag- 
netic aspect  relationship  for  different  geomagnetic  conditions  for  both  discrete 
and  diffuse  auroral  forms.  The  tapes  from  five  nights  of  data  collecting  during 
the  fall  of  1975  and  the  spring  of  1976  served  as  data  cases.  Specific  data  bases 
were  selected  with  the  aid  of  magnetometer  data  from  Churchill,  Canada,  as 
well  as  the  auroral  maps  and  histograms  which  had  been  generated  for  general 
support  of  the  various  auroral  study  tasks.  Measured  echoes  generally  had 
signal -to -noise  ratios  in  excess  of  15  dB  and  as  high  as  101  dB.  A noise  thresh- 
old required  detected  returns  to  exceed  about  10  dB  signal-to-noise  ratio,  thus 
limiting  errors. 


3.  1 The  Reflectivity  Scatter  Plot 

The  scatter  plot  was  chosen  as  the  media  for  outputting  the  processed 
aspect  data.  An  example  of  a typical  scatter  plot  is  shown  in  Figure  3-1.  These 
scatter  plots  are  generated  by  calculating  the  reflectivity  and  off -perpendicular 
angles  of  each  echo  value.  In  this  figure,  a linear  fit  to  the  data  has  been  per- 
formed manually. 

Two  important  factors  which  must  be  defined  are  the  calculation  of  re- 
flectivity and  the  off -perpendicular  angle.  A detailed  development  of  the  equation 
associated  with  these  parameters  has  been  provided  in  a previous  PAR  auroral 
study  report. [2],  The  definition  of  reflectivity  is  simply  the  effective  radar 
cross-section  per  unit  of  illuminated  volume.  This  is  determined  by  dividing 
the  apparent  radar  cross-section  <r  by  the  volume  of  the  range  resolution  cell 
V(R,p),  where  R is  the  apparen*  range  of  the  echo  and  p is  the  total  off-boresight 
angle.  The  PAR  is  a phased  array  radar,  the  beamwidth,  gain,  etc.  , depend 
on  the  off-boresight  angle.  These  dynamic  parameters  must  be  taken  into 
account  in  the  processing  of  each  echo. 


The  reflectivity  (Z)  is  given  by: 


V(R, P) 


2 / 3 

m / m 


From  [2],  the  effective  radar  cross-section  is  given  by: 


3 4 

(4r)  pR  R* 

2 3 

Prj>G>pGT^X  Lj^LpL<j>Cos  p 
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where: 


= the  echo  power  measured  at  the  signal  detector. 

P-p  = the  peak  power  measured  at  some  convenient  point  in 
the  system. 

L>P  = the  loss  from  the  transmit  power  measurement  point  to 
the  antenna  face. 

Grp  - the  antenna  power  gain  for  transmit. 

Gj^  = the  antenna  power  gain  for  receive. 

\ = the  operating  wavelength. 

R = the  range  to  the  center  of  the  range  resolution  cell. 

p = the  angle  between  the  antenna  bore  sight  and  the  direction 

of  the  range  resolution  cell. 


LR  = the  loss  between  the  antenna  face  and  the  signal  detector. 


the  two-way  attenuation  along  the  propagation  path. 


After  combining  the  constant  terms  and  assuming  a nominal,  two 
pagation  path  loss  (Lp)  of  0.4  dB,  the  equation  becomes: 


-way  pro- 


pr  R- 


Q Cos  p 


m2/m3 


3.3 


Calculation  of  the  volume  of  the  range  resolution  cell  depends  on  whether 
the  auroral  arc  fills  the  radar  beam.  If  the  arc  only  partially  fills  the  beam, 
the  cell  volume  is  given  by: 

Veffp  - R AH  !LJ_  lei<£) 

2 Cos  0 

where: 

A H = the  assumes  thickness  of  the 
auroral  arc. 


m3  3.4 


backscattering  portion  of  the 
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c = the  nominal  velocity  of  light. 

T = the  compressed  radar  pulse  width. 

B@(0°)  = the  azimuth  beamwidth  broadside  to  the  antenna. 

0 = the  steered  azimuth. 

The  cos  0 term  adjusts  for  the  azimuth  component  of  beam  broadening 
as  it  is  steered  off-broadside. 


Combining  these  equations  and  inserting  the  necessary  parameters  into 
the  reflectivity  equations  yields: 


R- 


"partial 


1.  18xl016  A H Cos  0 Cos2p 


2/  3 

m /m 


Reflectivity  is  normally  represented  in  logarithmic  form  with  the  reference 
being  1 M vM  . The  logarithmic  form  is  given  by: 


ZdB  = 10  Log  Z 


In  this  form,  the  quantity  Z dB  is  called  ZdBm2/m^  or  ZdBm_j. 

Generation  of  the  scatter  plot  also  requires  the  calculation  of  the  off-per- 
pendicular angle  for  each  data  point.  This  process  requires  a model  of  the 
earth's  magnetic  field.  The  model  used  in  this  study  was  developed  using  the 
"International  Geomagnetic  Reference  Field  1965.0."  [9]  The  implementation 
employs  a sixth-order  spherical  harmonic  solution;  the  48  terms  used  were 
updated  to  1975  validity  using  the  provided  rate  coefficients. 

3.2  Special  Techniques  Employed  to  Enhance  the  Aspect  Data 


A number  of  problems  and  potential  problems  were  encountered  during 
the  aspect  study.  These  problems  were  minimized  or  eliminated  using  a number 
of  techniques  which  will  be  described  in  the  following  paragraphs. 
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3.2.  1 Pulse -to -Pulse  Echo  Decorrelation 

It  has  been  observed  that  the  pulse -to -pulse  echo  intensity,  from  a given 
point  in  space,  fluctuates  very  rapidly,  becoming  completely  decorrelated  in 
about  10  ms  [6],  The  observed  distribution  of  echo  intensities  approximates  the 
theoretical  Rayleigh  distribution  shown  in  Figure  3-2.  The  use  of  data  distri- 
buted in  this  fashion  will  produce  a scatter  plot  with  a significant  dispersion, 
thus  degrading  the  confidence  interval  of  the  curve  fit.  The  quality  of  the 
aspect  data  has  been  improved  through  the  point -by -point  averaging  of  several, 
usually  6 to  10,  volume  scans.  In  this  process,  called  " super  scanning,  " a set 
of  synthetic  range  bins  are  created  for  each  pointing  angle,  then  all  the  echo 
peaks  from  the  corresponding  range /bins  of  the  different  scans  are  averaged. 
This  result  is  used  as  an  estimate  of  the  true  mean  echo  intensity.  It  can  be 
shown  that  a significant  improvement  is  realized  from  this  operation.  Figure 
3 -3a  shows  the  reflectivity  of  all  the  echo  peaks  from  a single  pulse  plotted 
versus  their  off -perpendicular  angles.  Figure  3-3b  shows  the  results  when  nine 
pulses  from  the  same  pointing  angle,  but  from  nine  different  scans,  have  been 
averaged.  The  improvement  in  dispersion  is  evident  immediately.  For  this 
reason,  averaging  or  superscanning  is  used  in  the  processing  of  all  aspect 
sensitivity  data  used  in  this  study. 

3.  2.  2 The  Effects  of  Resolution 

The  primary  factors  which  tend  to  degrade  the  results  of  aspect  sensi- 
tivity analysis  are  the  association  of  an  echo  sample  with  the  wrong  off-perpen- 
dicular  angle,  the  allowing  of  other  natural  variations  in  auroral  reflectivity 
to  be  attributed  to  changes  in  off-perpendicular  angle,  and  the  use  of  a geomag- 
netic field  model  which  is  not  representative  of  the  field  which  existed  during 
the  data  collection. 


The  errors  which  result  from  limitations  in  radar  resolution  have  a 
significant  effect  on  the  outcome  of  the  analysis.  The  types  of  resolution  which 
are  of  concern  include  not  only  the  range,  angle,  and  beamwidth  resolution  of 
the  radar,  but  temporal  resolution  determined  by  the  time  required  to  perform 
the  volume  scans.  The  faster  these  volume  scans  can  be  performed,  the  less 
large-scale  restructuring  takes  place  during  the  scan. 


Range  resolution  is  determined  primarily  by  the  effective  pulse  width 
of  the  radar.  Poor  range  resolution  not  only  results  in  the  association  of  echo 
peaks  with  the  wrong  off -perpendicular  angles,  but  may  result  in  an  overall 
distortion  in  the  shape  of  the  aspect  response  characteristic.  This  phenomena 
has  been  reported  by  two  sources  [6]  [7],  Figure  3-4a  shows  an  aspect  sensi- 
tivity plot  at  the  PARL  from  a 400  ps  continuous  wave  pulse  and  Figure  3-4b 
shows  the  effect  of  using  a chirped  pulse  with  a compres^d  pulse  width  of  6 ps. 
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SCATTER  PLOTS  SHOWING  EFFECTS  OF  SUPERSCANNING 


Single  Scan 


OFF-PERPEUDICULAP 


Figure  3 -3a 
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SCATTER  PLOTS  SHOWING  EFFECTS  OF  SUPERSCANNING 

(Continued) 


Nine  Scans 
Superscanned 
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Figure  3-3b 


In  this  case,  the  slopes  of  the  two  aspect  functions  are  similar,  but  the  absolute 
reflectivity  is  about  10dBm2/m3  less  for  the  CW  pulse.  This  phenomena  occurs 
because  in  each  range  resolution  cell,  the  power  measured  in  the  echoes  is  a 
composite  of  the  reflections  from  all  regions  within  this  cell;  the  total  cell  volume 
is  then  used  to  calculate  the  reflectivity.  If  there  are  inconsistencies  in  the 
scattering  medium,  then,  for  longer  pulse  widths,  the  measured  echo  may  no 
longer  be  representative  of  the  reflectivity  of  a finer-grained  irregularity  struc- 
ture. Because  of  its  increased  resolution,  the  shorter  pulse  width  produces  a 
more  accurate  representation  of  the  local  reflectivity. 

Angular  resolution  is  determined  by  the  accuracy  involved  in  pointing 
the  center  of  the  radar  beam.  This  error  is  usually  small  in  comparison  to  the 
effects  of  beamwidth  and  divergence.  The  beamwidth  resolution  problem  occurs 
because  in  the  search  mode,  the  radar  associates  the  angular  origin  of  all  echoes 
with  the  center  of  the  beam.  The  scatterers  from  which  the  measured  energy 
was  reflected  are  not  necessarily  located  at  the  center  of  the  beam.  In  fact, 
when  the  reflectivity  of  the  scatterers  is  high,  echoes  which  originated  several 
degrees  from  the  center  of  the  beam  can  still  be  detected,  as  shown  in  Figure 
3-5  which  shows  the  region  near  the  main  lobe  of  a radar  beam.  This  means 
that  during  the  radar's  volume  scan  the  same  scattering  region  may  be  detected 
several  times.  Detections  with  the  strongest  amplitude  will  occur  when  the  beam 
is  pointed  directly  at  the  scatterer,  but  weaker  echoes  will  occur  whenever  the 
beam  gets  within  a few  degrees  of  the  scattering  region.  This  effect  spreads 
and  distorts  the  apparent  aurora  in  a predictable  fashion.  For  the  PAR,  the 
maximum  off-beam-center  angle  at  which  the  radar  could  record  auroral  returns 
is  limited  by  the  sidelobe  blanking  system  or  Q -channel.  The  Q-channel  provides 
a high  probability  of  rejecting  most  echoes  which  originate  from  over  about  1.  7° 
from  beam  center.  It  should  be  noted  that  a number  of  sidelobe  returns  will 
also  be  detected  because  of  the  statistical  properties  of  the  auroral  echoes  and 
because  the  Q-channel  becomes  less  effective  for  weak  echoes. 

It  is  important  to  understand  that  the  echo  power  measured  at  any  instant 
by  the  receiver  is  a superposition  of  echo  power  from  the  same  slant  range  but 
from  all  directions  --  main  beam  as  well  as  sidelobes.  When  there  is  large- 
scale  blanketing  aurora,  significant  power  may  be  reflected  from  many  angles. 

The  Q-channel  cannot  eliminate  these  echoes  when  they  return  simultaneously 
with  echo  power  in  the  main  lobe.  This  unwanted  power  is  thus  added  to  and 
processed  with  the  main  beam  echoes.  An  unavoidable  error  in  the  determination 
of  auroral  reflectivity  has  thus  been  identified. 

The  divergence  of  the  conical  beam  makes  this  distortion  a function  of 
range,  thus  the  spreading  is  more  extreme  at  the  greatest  ranges.  Figure  3-6, 
a profile  plot  of  an  auroral  arc,  illustrates  this  effect  very  clearly.  The  appar- 
ent thickness  of  the  arc  is  a function  of  range  and  appears  to  be  converging  to 
its  true  value  directly  over  the  radar,  as  would  be  expected.  Unfortunately,  loss 
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Angle  Off  Array  Normal 


in  signal  strength  due  to  the  magnetic  aspect  response,  as  well  as  hardware 
limitations,  precluded  direct  overhead  observation.  But  with  some  analysis 
it  can  be  shown  that  the  arc  is  very  thin.  In  a region  of  high  reflectivity  such 
as  that  shown  in  Figure  3-6,  auroral  echoes  can  be  expected  to  originate  out  to 
the  extreme  of  the  1.  7°  blanking  limit.  The  weak  echoes  above  and  below  the 
dense  region  result  from  sidelobe  returns  which  are  not  blanked  by  the  Q- chan- 
nel. Thus,  near  the  radar  at  the  cutoff  range,  the  beamwidth  effect  should  in- 
crease the  apparent  thickness  of  the  arc  by  about  18  km.  The  "apparent"  thick- 
ness of  the  arc  is  shown  in  Figure  3-6  to  be  only  about  20  km.  It  is  apparent 
that  if  the  reflecting  region  was  much  over  a few  km  thick,  the  profile  would  be 
correspondingly  thicker  at  the  indicated  range.  From  this  simple  analysis,  the 
thickness  of  the  reflecting  region  could  be  estimated  to  be  from  2 to  5 km.  This 
estimate  is  not  far  from  the  results  of  previous  studies  [7,  10]  in  which  the  arc 
was  estimated  to  be  6 to  8 km  thick  in  the  first  case  and  in  the  latter,  2.  5 km 
thick.  Unger  determined  the  arc  thickness  by  deconvoluting  the  measured  data 
by  using  the  known  characteristics  of  the  radar  beam  [7]. 

Using  this  same  technique,  it  is  possible  to  calculate  the  shape  of  the 
distorted  arc  as  it  appears  to  the  radar.  The  boundary  of  this  region  has  been 
plotted  on  the  profile  plot  of  Figure  3-6.  It  can  be  seen  that  the  observed  arc 
matches  very  closely  the  predicted  shape. 

It  can  be  shown  by  this  analysis  of  beamwidth  effects  that  large  errors 
can  occur  if  one  attempts  to  use  the  ba^kscatter  data  from  all  altitudes  to 
characterize  the  aspect  response  function.  Figure  3-7,  a profile  view  of  the 
"apparent"  arc,  shows  the  off -perpendicular  contours  at  0°  azimuth.  It  can 
be  seen  from  this  plot  how  the  apparent  location  of  the  dispersed  echoes  changes 
their  relationship  to  the  magnetic  field.  This  results  in  an  improper  determin- 
ation of  the  aspect  response  function.  Scattering  centers  can  be  mislocated  by 
over  60  km  resulting  in  errors  of  over  2°  off -perpendicular.  It  is  therefore 
important  to  eliminate  any  echoes  which  did  not  originate  from  near  the  center 
of  the  beam.  This  can  be  done  easily  because  of  one  very  important  fact: 
when  an  echo  originates  near  beamcenter,  both  the  apparent  and  the  actual 
location  of  the  scatterer  conicide.  Furthermore,  if  the  altitude  of  the  arc  is 
known  and  only  echoes  which  originated  near  this  altitude  are  selected,  then 
essentially  all  off-axis  echoes  will  have  been  eliminated.  From  Figures  3-6 
and  3-7,  the  altitude  of  the  arc  can  be  closely  estimated  to  be  about  105  km. 

Another  means  of  estimating  the  altitude  of  the  center  of  the  arc  is  by 
use  of  a scatter  plot  to  generate  a profile  of  reflectivity  versus  altitude.  Figure 
3-8  has  been  prepared  to  illustrate  the  techniques  to  be  described.  The  radar 
beam  intersects  the  thin  auroral  arc  at  a low  elevation,  as  shown.  Because  of 
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PROFILE  VIEW  SHOWING  AURORAL  REFLECTIVITY 
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its  high  resolution,  there  are  many  range  resolution  cells  which  will  return 
echoes  to  the  radar.  The  radar,  however,  cannot  properly  locate  the  angle 
of  the  echoes  from  those  cells  which  are  not  on  beamcenter. 

The  radar  provides  an  accurate  measurement  of  range  to  the  cell.  How- 
ever, the  direction  of  the  cell  can  only  be  assumed  to  lie  along  the  central  angle 
of  the  beam  (the  dashed  line  in  Figure  3-8).  The  echo  cells  are  thus  mislocated 
in  altitude  except  for  the  cell  which  was  actually  in  the  center  of  the  beam. 
Furthermore,  the  gain  of  the  antenna  is  lower  for  all  but  the  center  cell,  so  the 
relative  response  during  a single  pulse  would  be  expected  to  resemble  that 
shown  in  the  lower  portion  of  Figure  3-8.  This  plot  peaks  near  the  location  of 
the  center  cell  which  is  also  the  cell  with  minimum  altitude  error  and  thus  truly 
represents  the  reflectivity  of  the  arc.  Thus,  on  a scatter  plot  showing  reflec- 
tivity versus  altitude,  a peak  should  occur  at  the  altitude  where  the  center  of 
the  arc  occurs.  Figure  3-9  is  a single  pulse  example  showing  this  phenomena 
and  that  the  arc  is  indeed  centered  near  105  km.  Figure  3-10  is  a similar 
plot  of  all  the  pulses  in  a scan.  The  shape  of  this  curve  demonstrates  the  dis- 
tortion in  the  "apparent"  reflectivity  which  occurs  if  echoes  are  used  which  did 
not  originate  near  beamcenter.  A plot  of  reflectivity  versus  off -perpendicular 
angle  for  a single  pointing  angle  is  shown  in  Figure  3-11.  Data  from  all  altitudes 
is  included  in  this  plot.  In  this  plot,  the  off-perpendicular  angle  increases  with 
slant  range  and  thus  with  increasing  altitude.  From  the  plot,  it  can  be  seen 
that  the  reflectivity  steadily  increases  with  altitude  until  it  peaks;  it  then  falls 
above  this  height.  It  can  be  shown  that  this  peak  occurs  near  the  center  of  the 
arc,  in  this  case  about  105  km.  The  shape  of  the  curve  approximates  the  shape 
of  the  peak  of  the  radar  beam,  in  that  the  image  of  the  aurora  is  the  convolution 
of  the  entire  radar  beam  with  the  "actual"  auroral  scattering  region.  The  curve 
trails  off  at  a decreased  slope  beyond  4°  off -perpendicular  because  of  the  super- 
position of  side  lobe  returns  with  those  from  the  mainbeam. 

Figure  3-12,  which  includes  a large  number  of  pointing  angles,  shows 
wrat  happens  if  auroral  backscatter  data  from  all  altitudes  is  used  in  the  scatter 
plot.  The  curves  from  the  individual  pointing  angles  are  still  discernable  as 
well  as  the  low  reflectivity  independent  side  lobe  returns. 

It  is  readily  apparent  from  this  plot  that  little  confidence  could  be  placed 
in  any  curve  fit  through  this  superposition  of  data.  Figure  3-13  has  been  gen- 
erated using  the  same  data  as  Figure  3-12,  except  that  data  from  105  to  106  km 
only  is  included.  It  can  be  seen  that  these  data  points  consist  of  the  peak  re- 
flectivities from  the  individual  pulses  or  from  near  the  center  of  the  reflecting 
arc.  Asa  result  of  this  analysis,  all- data  for  succeeding  scatter  plots  was 
selected  from  within  1 km  of  the  estimated  center  of  the  reflecting  region.  The 
peak  of  the  altitude  versus  reflectivity  plots  is  broad  enough  so  that  an  error 
of  even  +3.5  km  in  the  estimation  of  this  peak  would  result  in  little  error. 
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ALTITUDE  VERSUS  REFLECTIVITY  PROFILE  FOR  AN  ENTIRE  SCAN 
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Another  problem  encountered  in  the  aspect  study  results  from  regional 
variations  in  auroral  reflectivity.  These  inconsistencies  which  appear  as  voids, 
edges,  or  overall  shifts  in  reflectivity  are  not  caused  by  the  normal  aspect  re- 
lationship. An  example  of  an  auroral  scan  with  gaps  is  shown  in  Figure  3-14. 

The  off-perpendicular  contours  have  been  overlaid  on  this  data  providing  a 
straightforward  method  for  determining  what  magnetic  aspects  will  be  represented 
on  the  scatter  plot.  This  figure  show3  the  sharp  regional  changes  in  reflectivity 
which  are  often  observed.  The  shifts  in  the  overall  reflectivity  for  another  such 
scan  are  apparent  in  Figure  3-15.  The  scatter  plot  has  a pattern  containing 
three  segments,  each  of  which  is  obeying  the  same  predictable  aspect  response 
function  except  that  each  is  offset  by  an  overall  scale  factor.  Further  investi- 
gation showed  that  each  segment  of  this  scatter  plot  was  generated  from  echoes 
originating  in  widely  separated  regions. 

In  order  to  minimize  this  problem,  care  must  be  taken  to  select  data 
from  auroral  forms  which  are  as  homogeneous  as  possible.  This  selection  in- 
creases the  probability  that  the  changes  in  apparent  reflectivity  observed  on 
the  scatter  plot  truly  represent  the  effects  of  the  aspect  response  function.  As 
a result  of  these  observations,  each  scan  used  in  the  PAR  aspect  study  was  first 
plotted  on  a top-down  map  with  off-perpendicular  contour  overlays  as  was  shown 
earlier  in  Figure  3-14,  and  the  data  used  to  generate  the  scatter  plots  was  taken 
from  regions  in  which  the  auroral  echoes  appeared  homogeneous.  An  important 
factor  was  to  assure  that  the  data  was  selected  from  a region  including  the 
widest  range  of  off-perpendiculars  within  the  smallest  possible  region.  This 
minimized  the  degradation  in  the  results  due  to  the  regional  variations  in  auroral 
makeup  which  were  discussed  previously.  This  data  selection  was  accomplished 
with  the  aid  of  a specially  designed  spatial  filter  and  a special  top-down  overlay 
which  shows  the  position  at  which  each  beam  penetrates  the  altitude  shell  at  the 
center  of  the  arc.  Using  this  overlay,  shown  in  Figure  3-16,  it  was  possible 
to  select  the  exact  steering  angles  from  which  to  obtain  data.  By  observing  the 
top-down  slice  corresponding  to  the  105  km  altitude,  it  was  possible  to  pick 
the  beam  positions  from  which  echoes  were  detected  at  extreme  off-perpendicular 
angles  while  assuring  elimination  of  voids,  edges,  etc. 

One  more  problem  inherent  in  this  or  any  aspect  study  is  determining 
whether  the  magnetic  field  model  realistically  represents  the  actual  magnetic 
conditions  which  were  valid  at  the  time  of  that  collection.  Associated  with  an 
auroral  arc  is  a strong  current  flow  which  distorts  the  local  geomagnetic  field. 

Two  techniques  were  employed  to  minimize  field  distortion- related 
errors  for  discrete  arcs.  The  first  was  to  examine  the  magnetometer  data  from 
Fort  Churchill,  Canada,  participating  research  station  located  near  the  auroral 
arc  being  observed.  The  location  of  Fort  Churchill  is  shown  on  the  map  in 
Figure  1-1.  Backscatter  data  was  not  used  for  aspect  calculations  when  the 
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geomagnetic  field  was  unstable  or  if  any  of  the  components  had  deviated  by  much 
more  than  100  gammas  from  the  nominal  baseline. 


The  second  technique  is  based  on  the  fact  that  the  auroral  arc  is  normally 
aligned  with  the  magnetic  L-shells.  This  phenomenon  can  be  used  as  an  indicator 
of  the  magnetic  distortion  by  plotting  the  L-shells  on  a top-down  radar  map  of 
the  arc  as  shown  in  Figure  3-17.  In  this  study,  only  well-aligned  arcs  have  been 
used  to  calculate  the  aspect  response.  Although  this  technique  alone  does  not 
assure  a minimization  of  effects  from  magnetic  distortion,  combining  it  with 
the  magnetometer  screening  offers  a measure  of  improvement  in  results,  espe- 
cially in  repeatability  over  the  use  of  randomly  selected  data.  The  results  of 
the  analysis  reveal  a phenomena  associated  with  discrete  arcs  which  suggest 
that,  although  the  above  technique  may  assure  a degree  of  repeatability,  a signi- 
ficant local  magnetic  distortion  is  inherent  with  these  discrete  forms.  It  appears 
that  this  distortion  may  result  in  an  extreme  shift  in  the  magnetic  aspect  response 
characteristic.  This  phenomena  will  be  discussed  further  in  a later  section. 
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TOP-DOWN  REFLECTIVITY  MA  P WITH  L-SHELL  CONTOUR  AT  105  km 


4. 


DATA  CASES  USED  IN  THE  ANALYSIS 


Several  test  cases  were  carefully  selected  for  the  aspect  study.  They 
consist  of  both  diffuse  and  discrete  forms.  A special  distinction  has  been  made 
between  these  two  basic  forms  because  of  the  recognition  of  a clear  difference 
between  the  ’'apparent"  aspect  response  characteristics  of  each  type.  Data  rep- 
resented in  this  study  was  collected  on  September  27,  1975,  and  on  March  10, 

18,  26,  and  31,  1976.  The  next  few  sections  show  the  results  of  the  data  reduc- 
tion for  each  test  case.  Included  is  a top-down  auroral  map  with  overlaid  105  km 
off-peipendicular  contours  showing  the  region  where  the  scatter  olot  data  wf  s 
selected.  Two  types  of  scatter  plots  are  generally  included.  The  first  shows  a 
reflectivity  versus  altitude  plot  used  to  locate  the  center  of  the  reflecting  arc 
as  described  in  Section  3;  the  second  shows  the  reflectivity  vf  sus  off-perpendic - 
ular  angle.  A curve -fit  has  been  performed  on  each  of  these  scatter  plots. 

4.  1 Test  Case  270/1058 

The  first  test  case  occurred  on  September  27,  1976.  A top-down  auroral 
reflectivity  map  of  this  scan  is  shown  in  Figure  4-1.  This  auroral  scan  occurred 
during  a period  of  relatively  low  magnetic  activity  as  indicated  by  the  nearby 
Fort  Churchill,  Canada,  magnetometer  data.  The  deviations,  from  their  cali- 
brated baseline  of  the  three  magnetic  components,  were  80  ? for  the  X component, 
45?  for  the  Y component,  and  5?  for  the  Z,  component.  The  three  traces  were 
essentially  flat  during  the  time  of  the  scan.  The  riometer  reading  was  also  quiet 
and  at  a null.  The  peak  signal-to-noise  ratio  of  these  echoes  was  60  dB  and  they 
averaged  about  30  dB.  The  map  shows  a typical  quiet  diffuse  form  or  possibly 
the  south  edge  of  an  arc.  Good  aspect  data  may  be  expected  due  to  low  magnetic 
field  distortion  within  the  region  covered.  The  reflectivity  versus  altitude  plot 
in  Figure  4-2  shows  the  center  of  the  arc  to  be  near  105  km,  so  aspect  sensi- 
tivity data  was  selected  only  from  within  the  105-106  km  band.  By  examining  the 
relationship  of  the  off -perpendicular  contours  to  the  aurora,  shown  in  Figure  4-1, 
it  can  be  seen  that  the  off -perpendicular  angles  from  which  echoes  originated 
range  from  2°  to  8°.  To  determine  the  aspect  sensitivity.  Figure  4-3  was 
generated  using  echo  data  from  the  small  region  A shown  in  Figure  4-1.  Data 
for  the  second  scatter  plot  (Figure  4-4)  came  from  region  B.  These  two  cuts 
were  made  in  order  to  evaluate  any  differences  in  the  apparent  aspect  response 
when  data  was  sampled  approximately  parallel  to  the  L-shell  as  in  region  A or 
perpendicular  to  the  L-shells  as  in  region  B. 

The  two  scatter  plots  demonstrate  a clear  difference  in  apparent  aspect 
sensitivity.  When  the  data  points  are  selected  along  region  A,  a 7.  34  dB/deg 
slope  was  measured,  but  from  the  region  perpendicular  to  the  L-shells  (region 
B),  a 16.3  dB/deg  slope  was  found. 
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Test  Case  70/666 


This  test  case  involves  a discrete  arc  as  shown  in  Figure  4-5.  This  arc 
consists  of  the  main  reflecting  band  shown  in  the  lower  portions  of  the  figure 
and  a smaller  band  above  it.  The  only  magnetic  data  available  was  at  14:14  UT 
when  the  Ap  index  was  about  33  and  a minor  geomagnetic  storm  was  in  progress. 
The  absence  of  magnetometer  data  from  Fort  Churchill  makes  it  impossible  to 
estimate  the  instantaneous  geomagnetic  distortion;  however,  there  is  good  align- 
ment of  the  arc  with  the  L-shell  contours.  The  peak  s/n  was  54  dF,  the  average 
about  24  dB.  Two  scatter  plots  were  generated  using  data  from  the  105-106  km 
band  as  indicated  by  the  peak  of  the  altitude  versus  reflectivity  profile  of  Figure 
4-6.  The  scatter  plot  shown  in  Figure  4-7  was  generated  from  data  selected 
from  region  A in  Figure  4-5.  The  scatter  plot  shown  in  Figure  4-8  resulted  from 
data  selected  in  region  B.  A distinct  difference  in  the  aspect  response  charac- 
teristics of  these  two  plots  is  readily  apparent.  One  unusual  observation  is  the 
very  high  (20  dB/deg)  slope  at  the  low  off-perpendicular  angles  as  seen  in 
Figure  4-7. 

4.  3 Test  Case  78/12 


This  is  a homogeneous  and  clearly  defined  discrete  arc  from  which  echoes 
were  measured  with  signal-to-noise  ratios  as  great  as  64  dB  and  averaging  about 
46  dB.  The  magnetic  activity  was  low  with  an  ap  of  only  16  reported  at  14:14  UT. 
Figure  4-9  shows  that  echoes  were  detected  throughout  the  range  of  angles  from 
2°  to  8°  off-perpendicular.  Figure  4-10  shows  the  center  of  the  arc  at  about 
105  km.  The  aspect  response  curve  shown  in  Figure  4-11  was  created  from 
data  parallel  to  the  L-shells  (region  A in  Figure  4-9).  The  curve  in  Figure  4-12 
was  created  from  echoes  along  region  B. 

The  slopes  of  these  curves,  which  were  manually  fit  to  the  data,  reveal 
again  the  gross  differences  noted  when  attempting  to  determine  aspect  sensitivity 
from  different  sets  of  data.  The  differences  in  the  shape  of  these  aspect  response 
curves  will  be  discussed  in  more  detail  in  a later  section,  but  the  differences  are 
felt  only  to  be  apparent,  but  are  in  fact  the  result  of  magnetic  distortion  from 
the  current  flow  in  the  discrete  arc. 
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Figure 


4.  4 Test  Case  86/576 

The  scan  representing  this  test  case  is  shown  in  Figure  4-13.  The  map 
shows  this  discrete  arc  including  off -perpendicular  angles  from  2°  to  8°.  The 
Churchill  magnetometer  was  steady,  showing  a deviation  from  baseline  of  +60  + 
for  the  X component,  +5?  for  the  Y component,  and  -80+  for  the  Z component. 
The  arc  was  well-aligned  with  the  L-shells  as  shown  in  Figure  4-14.  Figure 
4-15  shows  the  arc  to  be  centered  near  105  km. 

Sample  regions  parallel  and  perpendicular,  respectively,  to  the  L-shell, 
shown  as  region  A and  region  B in  Figure  4-13,  yield  the  scatter  plots  given  in 
Figures  4-16  and  4-17.  The  contrast  between  the  apparent  aspect  sensitivity 
of  the  two  samples  is  again  evident.  A slope  of  only  5.83  dB/deg  was  calculated 
for  the  region  approximately  parallel  to  the  L-shell  contour  (region  A)  and 
11.4  dB/deg  for  the  data  from  the  perpendicular  region  (region  B).  The  curve- 
fit  used  in  the  figure  required  manual  screening  of  the  points  on  either  end  of 
the  scatter  plot  because  the  severe  reflectivity  dropoff  noted  on  either  end  is  a 
direct  result  of  the  sample  region  crossing  the  arc  boundaries  shown  as  C and 
D on  Figure  4-13.  At  these  boundaries  there  is  a sharp  transition  from  rela- 
tively strong  reflectivity  to  zero  reflectivity.  This  results  in  the  effects  noted 
and  is  especially  clear  near  the  2°  off -perpendicular  angle  in  Figure  4-13. 

These  edge  effects  are  not  representative  of  the  aspect  response  and  must  be 
manually  detected  and  eliminated  in  the  curve -fit  and  analysis  process. 

4.  5 Test  Case  86/632 

This  discrete  arc,  shown  in  Figure  4-18,  occurred  during  quiet  magne- 
tic conditions  nearly  identical  to  the  previous  test  case.  The  two  scatter  plots 
were  generated  from  the  two  cuts  along  region  A and  region  B and  from  the 
105-106  km  altitude  band  as  indicated  by  the  altitude  profile  of  Figure  4-19. 

The  scatter  plot  in  Figure  4-20  was  generated  from  data  within  region  A of 
Figure  4-18,  or  generally  parallel  to  the  L-shells.  The  scatter  plot  of  Figure 
4-21  was  generated  from  data  along  region  B,  perpendicular  to  the  L-shells. 

It  should  be  noted  that  only  three  scans  were  available  for  averaging. 

This  fact  is  evident  in  the  more  than  usual  dispersion  of  the  scatter  plots, 
especially  notable  in  Figure  4-20.  The  curve -fit  to  the  scatter  plot  of  Figure 
4-21  required  manual  correction  to  eliminate  the  edge  effects.  By  studying 
how  the  off -perpendicular  contours  fit  the  map  on  Figure  4-18,  it  can  be  seen 
that  only  the  region  from  about  2.5°  to  3.  75°  off -perpendicular  can  yield  echoes 
free  of  edge  effects. 

Again,  the  sample  region  B perpendicular  to  the  L-shell  produces  a 
steeper  aspect  response  slope  than  does  the  data  from  region  A. 
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MAP  OF  SCAN  576  SHOWING  L-SHELL  CONTOURS  AT  105  KM 


Reproduced  from 
best  available  copy. 
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4.6 


Test  Case  86/5557 


As  indicated  in  Figure  4-22,  this  map  shows  a large  diffuse  auroral  form. 
The  poleward  limits  of  the  reflectiong  region  are  determined  by  the  radar  hori- 
zon, and  the  boundary  near  the  radar  results  from  the  minimum  range  of  the 
sensor.  Reports  from  the  Finley  Air  Force  Radar  Station,  south  of  the  PAR, 
confirmed  that  the  equatorial  edge  of  the  arc  extended  far  south  of  the  PAR. 
Magnetic  activity  at  Fort  Churchill  was  very  steady  at  this  time.  The  magne- 
tometer was  steady  and  showed  a deviation  of  -407  in  the  X component,  -75 7 
in  the  Y component,  and  -1407  in  the  Z component. 

The  center  of  the  reflectivity  arc  is  located  near  105  km  as  shown  by 
Figure  4-23.  The  altitude  profile  is,  however,  unusual  because  of  the  gaps  on 
either  side  of  the  primary  plot.  No  attempt  will  be  made  to  explain  this  phe- 
nomenon. Three  aspect  response  scatter  plots  have  been  generated,  one  from 
each  of  the  three  regions,  labeled  A,  B,  and  C in  Figure  4-22.  The  first 
scatter  plot.  Figure  4-24,  comes  from  region  A;  the  slope  is  5.5  dB/deg. 

The  second  scatter  plot,  shown  in  Figure  4-25,  was  generated  with  data  from 
region  B;  its  slope  is  4.  84  dB/deg.  The  final  scatter  plot  was  generated  from 
a combination  of  data  from  regions  B and  C.  Region  C includes  data  from  the 
maximum  attainable  off -perpendicular  angles.  The  resultant  scatter  plot  of 
Figure  4-26  covers  a wide  range  of  aspect  angles  and  shows  how  the  slope 
flattens  at  high  aspect  angles.  This  effect  was  described  earlier  and  has  been 
noted  by  others  [6,9]. 

This  test  case  is  the  first  clearly  diffuse  large-scale  aurora  investigated 
in  the  study.  It  should  be  noted  that  for  the  first  time  the  aspect  response  slope 
was  lower  when  the  data  was  selected  perpendicular  to  the  L-shells,  than  when 
selected  parallel  to  them,  as  evidenced  by  Figures  4-24  and  4-25.  In  fact,  there 
are  no  significant  differences  in  the  two  curves. 

4.  7 Test  Case  86/5571 


The  data  set  which  generated  the  auroral  map  in  Figure  4-27  was  collected 
near  in  time  to  the  previous  scan  and  the  apparent  geomagnetic  conditions  were 
nearly  identical.  The  arc  was  assumed  to  be  centered  near  105  km.  The  three 
scatter  plots,  shown  in  Figures  4-28,  4-29,  and  4-30,  were  generated  with  data 
from  regions  A,  B,  and  C,  respectively  (see  Figure  4-27).  Comparing  Figures 
4-28  and  4-29  again  shows  the  aurora  appears  more  aspect- sensitive  when 
sampling  perpendicular  to  the  L-shells  than  when  selecting  data  along  a line 
generally  parallel  to  the  contours.  The  scatter  plot  in  Figure  4-30  was  generated 
from  region  C data  and  displays  very  little  aspect  sensitivity.  All  the  curves 
demonstrate  a decreasing  slope  with  increasing  aspect  angle. 
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4.  8 Test  Case  86/5622 


This  scan  shows  a large  blanketing  diffuse  form,  very  much  like  the  pre- 
vious examples.  The  index  was  234  during  the  observation,  but  the  magne- 
tometer at  Fort  Churchill  was  steady  showing  about  the  same  deflections  as  test 
case  86/555  7.  The  off-perpendicular  contours  in  Figure  4-31  show  again  echoes 
detected  from  over  a wide  range  of  magnetic  aspect  angles,  out  to  about  18°. 
Figure  4-32  shows  the  center  of  the  arc  to  be  located  near  the  center  of  the 
105-106  km  altitude  band,  from  which  all  data  for  the  three  scatter  plots  was 
selected. 

The  first  scatter  plot,  shown  in  Figure  4-33,  was  created  from  data 
taken  along  region  A in  Figure  4-31.  Due  to  its  distorted  appearance,  the  slope 
is  given  at  only  one  point.  The  scatter  plot  shown  in  Figure  4-34  was  created 
with  the  data  from  region  B,  and  the  third  scatter  plot,  shown  in  Figure  4-35, 
represents  the  aspect  response  obtained  by  combining  data  from  regions  B and 
C.  This  aspect  response  demonstrates  a steady  decrease  in  slope  from  about 
6 dB/deg  to  1.9  dB/deg.  Although  the  last  two  scatter  plots  demonstrate  a 
well-behaved  aspect  response,  there  was  apparently  a high  level  of  magnetic 
field  distortion  aligned  in  some  way  with  the  L-shells  which  resulted  in  the 
distortion  seen  in  Figure  4-33.  Explaining  this  phenomena  is  beyond  the  scope 
of  this  report,  but  such  an  effect  could  result  from  a localized  current  com- 
ponent directed  perpendicular  to  the  L-shell  contours. 
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4.  9 Test  Case  92/681 

Although  diffuse,  the  structure  of  the  reflecting  region  was  changing 
quite  rapidly  during  the  10  volume  scans  used  in  this  test  case.  This  fact  is 
evident  from  the  scans  in  Figures  4-36  and  4-37,  the  scans  at  the  beginning  and 
end  of  the  observation  period.  Unfortunately,  magnetometer  data  from  Fort 
Churchill  was  not  available,  and  the  ap  index  was  over  200,  so  the  local  magne- 
tic conditions  are  not  clear.  However,  at  the  PAR,  diffuse  aurora  has  usually 
been  noted  only  during  steady  local  magnetic  activity. 

Figure  4-38  shows  the  reflecting  arc  to  be  centered  near  105  km  as 
indicated  by  the  altitude  versus  reflectivity  plot.  The  first  aspect  response 
curve  of  Figure  4-39,  created  from  data  in  region  A,  is  given  in  Figure  4-36. 
This  sample  region  is  generally  parallel  to  the  magnetic  L-shells  and  demon- 
strates almost  no  aspect  sensitivity  up  to  about  5°  off -perpendicular  where  the 
slope  rapidly  changes  to  nearly  12  dB/deg.  More  evidence  of  this  unusual  phe- 
nomena can  be  seen  in  the  top-down  map  of  Figure  4-37.  Notice  the  strongly 
reflecting  central  region  and  the  sudden  dropoff  in  reflectivity  starting  between 
4°  and  5°  off-perpendicular.  It  was  immediately  suspected  that  this  may  have 
resulted  from  saturation  of  the  radar  receiver,  but  subsequent  checks  assured 
that  no  suc'a  problem  occurred.  The  apparent  aspect  response  was  thus  properly 
represented.  Of  course,  this  does  not  eliminate  the  effects  of  magnetic  field 
distortion  during  this  extremely  active  auroral  occurrence. 

The  second  aspect  response  curve  is  associated  with  region  B,  perpen- 
dicular to  the  L-shell.  This  plot,  shown  in  Figure  4-40,  demonstrates  a more 
conventional  characteristic,  but  it  still  shows  an  apparent  increase  in  aspect 
sensitivity  at  greater  off-perpendicular  angles.  The  third  scatter  plot  includes 
data  from  regious  B and  C.  The  response  curve,  shown  in  Figure  4-41,  has 
a well-behaved  slope  out  to  about  12°  off-perpendicular  where  it  flattens  and 
even  appears  to  become  positive.  This  flattening  or  slope  reversal  has  been 
observed  in  other  test  cases  and  may  be  a result  of  radar  limitations  rather 
than  an  actual  change  in  aspect  sensitivity. 


The  final  aspec*  response  curve,  shown  in  Figure  4-42,  was  generated 
from  data  in  regions  B and  D.  Aspect  data  is  shown  covering  angles  from  2°  to 
about  17°  off-perpendicular.  The  slope  shift  at  about  5°  is  still  noticeable, 
and  the  data  distribution  becomes  less  well-defined  at  angles  from  7°  to  17°. 
Most  of  these  data  points  came  from  region  D which  was  apparently  somewhat 
perturbed  at  that  time. 
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4.  10  Test  Case  92/723 


This  test  case  has  been  included  because  the  scan  includes  echoes  with 
signal-to -noise  ratios  near  100  dB.  These  may  be  among  the  strongest  auroral 
echoes  ever  detected  by  a radar.  The  strongest  echoes  measured  by  PAR  some- 
what earlier  exceeded  101  dB  signal-to -noise  ratio  with  a reflectivity  of  about 
-TOdBj^/^S.  The  distribution  of  reflectivities,  after  11  scans  were  averaged, 
is  shown  in  the  histogram  in  Figure  4-44.  This  histogram  was  generated  with 
data  from  region  B in  Figure  4-43. 

The  ap  index  was  very  high  at  this  time  (#200)  and  no  local  geomagnetic 
data  is  available.  The  scatter  plots  for  this  test  case  were  generated  from  the 
105-106  km  altitude  band.  The  first  of  these  was  created  with  data  from  along 
region  A of  Figure  4-43.  The  data  from  this  region,  which  is  generally  parallel 
to  the  L-shells,  results  in  the  scatter  plot  in  Figure  4-45.  Evidence  of  some 
type  of  distortion  is  apparent  at  about  2°  off-perpendicular,  but  the  function  is 
quite  well-defined  at  other  aspect  angles  and  has  a slope  of  6.  68  dB/deg.  The 
second  aspect  response  function,  shown  in  Figure  4-46,  was  generated  with  data 
from  regions  B and  C.  More  distortion  is  evident  in  this  data,  but  the  slope  of 
the  aspect  response  was  clearly  quite  low  and  measured  about  2.  25  dB/deg.  One 
point  of  special  interest  in  this  case  is  that  the  apparent  slope  of  the  aspect  did 
not  decrease  as  drastically  as  in  other  cases.  This  may  be  because  the  reflec- 
tivity was  greater  than  for  other  cases  so  that  there  was  a high  mainbeam-to- 
sidelobe  power  ratio.  This  supports  the  author's  theory  that  the  apparent  de- 
crease in  slope  may  be  due  to  sidelobe  returns  as  will  be  discussed  in  Section  5. 
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Even  though  very  little  discussion  was  included  with  the  test  cases, 
many  of  the  indications  are  obvious  upon  observation  of  the  data.  The  most 
salient  point,  however,  is  that  in  spite  of  all  the  safeguards  employed  to  mini- 
mize the  ambiguities  in  the  processing  of  auroral  magnetic  aspect  data,  an 
array  of  aspect  response  curves  has  been  generated.  The  primary  informa- 
tion from  the  test  cases  has  been  compiled  in  Table  5-1,  which  clearly  shows 
some  of  the  effects  noted.  It  is  evident  from  the  individual  test  cases  that  the 
mechanism  responsible  for  the  apparent  aspect  sensitivity  relationship  is  a 
highly  dynamic  phenomena.  The  marked  differences  in  the  aspect  response 
curves  from  the  different  sample  regions  suggest  that  there  may  be  only  a very 
small  volume  of  validity  for  the  true  aspect  response  characteristics.  Depending 
on  the  sampling  method,  each  scan  will  probably  yield  innumerable  aspect  curves. 

It  is  clear  that  discrete  arc  forms  "appear"  to  have  a steeper  aspect 
sensitivity  slope  than  do  diffuse  forms.  This  phenomena  may,  however,  be  only 
an  apparent  effect,  since  the  current  systems  associated  with  discrete  arcs 
can  create  distortion  in  the  local  magnetic  field,  which  is  somewhat  different 
than  that  for  diffuse  aurora.  It  appears  that  the  mechanism  which  results  in 
aspect  sensitivity  may  be  more  complex  than  originally  thought  because  of  the 
apparent  difference  in  the  aspect  sensitivity  relationship  of  samples  from  differ- 
ent regions  of  a scan.  Another  reason  for  suspecting  a more  complex  aspect 
phenomena  is  that  in  a number  of  cases,  the  magnetic  field  distortion  appeared 
very  low,  and  calculations  based  on  available  magnetometer  data  predict  that 
the  changes  in  apparent  aspect  sensitivity  should  be  much  smaller  than  those 
observed.  This  effort  has  been  most  marked  when  comparing  aspect  response 
plots  generated  in  one  case  with  data  selected  parallel  to  the  L- shells  contours 
and  in  the  other  case,  perpendicular  to  these  contours. 

One  very  interesting  observation  is  that  in  every  observed  case,  the 
peak  reflectivity  occurs  at  almost  precisely  the  same  altitude.  This  phenomena 
would  seem  most  unexpected  in  nature. 

It  has  been  noted  ir*  most  applicable  test  cases  that  the  slope  of  the 
aspect  response  function  appears  to  flatten  out  at  large  off -perpendicular  angles.  i 

This  observation  must  be  viewed  with  caution,  for  it  may  be  simply  an  effect 
caused  by  the  limitations  of  the  radar.  The  radar  beam  consists  not  only  of 
the  main  beam,  but  numerous  sidelobes  pointing  in  many  directions.  The 
energy  measured  in  each  echo  is  the  superposition  of  echoes  from  the  scatterers 
in  the  main  beam  as  well  as  in  the  directions  of  all  sidelobes.  The  radar  asso- 
ciates all  of  this  echo  power  with  the  main  beam  echo.  The  echo  strength  is 
thus  increased  somewhat  by  energy  received  in  these  unwanted  echoes.  Even 
when  the  main  beam  power  is  low,  significant  side  lobe  power  may  still  be  detected 
and  processed  as  an  "apparent"  main  beam  echo.  The  sidelobe  blanking  system 
cannot  eliminate,  reduce,  or  affect  in  any  way  sidelobe  echoes  when  they  are 
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superimposed  on  detectable  main  beam  echoes.  This  problem  is  especially 
severe  in  large-scale  diffuse  aurora  because  there  is  reflecting  aurora  in  all 
directions,  therefore,  many  sidelobes  will  be  receiving  and  contributing  un- 
wanted power  to  the  main  beam  echoes. 

These  sidelobe  echoes  affect  the  measurement  of  the  aspect  response 
function  in  a predictable  fashion.  Where  the  main  beam  makes  small  off -per- 
pendicular angles,  the  total  superimposed  sidelobe  power  is  small  in  relation 
to  the  main  beam  power;  thus,  the  associated  distortion  to  the  aspect  curve  is 
small.  As  the  main  beam  moves  to  greater  off -perpendicular  angles,  where  it 
is  attenuated  by  aspect  sensitivity,  a number  of  sidelobes  will  always  be  pointed 
in  the  direction  of  small  aspect  angles.  The  echo  power  in  this  sidelobe  is  thus 
enhanced  and  the  ratio  of  side  lobe -to -main  beam  power  increases.  The  primary 
effect  is  that  the  aspect  response  slope  flattens  out,  and  at  large  aspect  angles, 
where  most  of  the  signal  power  is  actually  contributed  from  the  sidelobes,  the 
scatter  plot  approaches  a constant  mean  value.  In  one  test  case,  92/723,  the 
aspect  response  had  still  not  completely  flattened  even  at  17°  off-perpendicular. 
This  may  have  happened  because  the  main  beam  auroral  reflectivity  was  so  high 
at  the  large  aspect  angles  that  there  was  still  enough  power  from  the  main  beam 
to  offset  the  unwanted  sidelobe  power. 

Returning  to  the  analysis  of  the  aspect  response,  it  should  be  explained 
that  no  attempt  has  been  made  to  generate  an  average  or  composite  aspect 
response  curve  for  these  cases  because  of  the  widely  varying  results  already 
noted.  Much  has  already  been  done  to  assure  that  the  individual  aspect  curves 
represent  the  most  accurate  results  possible  using  the  existing  quiescent  magne- 
tic field  model.  It  is  apparent  that  some  changes  are  taking  place  in  either  the 
magnetic  field  or  the  auroral  reflection  mechanism  which  is  causing  these  aspect 
curves  to  vary.  Averaging  of  these  curves  would  only  mask  the  actual  phenomena. 
It  would  appear  that  much  more  could  be  gained  by  utilizing  the  difference  in 
these  curves  to  help  determine  more  accurately  the  mechanism  and  dynamics  of 
this  wave  action  which  makes  auroral  backscatter  possible.  Other  authors  [3,6] 
have  noted  that  aspect  sensitivity  varies  according  to  the  magnetic  index,  as  was 
shown  in  Figure  2-3.  It  has  been  observed  that  the  auroral  forms  and,  of  course, 
the  level  of  magnetic  distortion  also  varies  with  the  magnetic  indices.  For  this 
reason,  it  cannot  be  assumed  that  the  true  aspect  sensitivity  relation  varies. 

It  may  be  simply  that  the  quiescent  magnetic  field  model  may  no  longer  be  rep- 
resentative. These  indices  may  then  merely  provide  an  indicator  of  the  validity 
of  the  field  model  used. 

Another  of  the  differences  between  diffuse  and  discrete  aurora  can  be 
seen  by  comparing  the  differences  in  apparent  aspect  sensitivity  when  sampling 
the  two  generally  orthogonal  regions  where  the  sample  region  was  elongated 
paraLlel  to  the  L-shells  and  as  in  the  other  case,  perpendicular  to  them.  It  can 
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be  readily  ascertained  from  the  data  which  has  been  given  that,  for  these 
examples,  the  apparent  aspect  sensitivity  was  always  greater  for  discrete 
aurora  when  the  data  was  sampled  from  a region  elongated  perpendicular  to  the 
L-shell.  For  diffuse  aurora,  aspect  sensitivity  was  normally  greater  when  the 
data  was  sampled  from  a region  generally  parallel  to  the  L-shells.  The  effect 
was,  however,  less  marked  than  for  discrete  aurora.  The  phenomena  just 
described  is  not  too  surprising  when  one  considers  the  differences  in  the  magne- 
tic field  components  which  are  induced  by  the  current  system  associated  with 
the  two  types  of  auroral  forms.  A further  discussion  of  these  induced  fields 
is  included  in  Section  6. 
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6. 


PERTURBATION  OF  THE  GEOMAGNETIC  FIELD  BY  IONOSPHERIC 
CURRENTS 


It  has  been  observed  that  there  is  considerable  variation  in  the  apparent 
aspect  sensitivity  of  the  different  test  cases,  particularly  between  diffuse  and 
discrete  forms.  Also,  aspect  sensitivity  measured  parallel  to  the  L-shell  differs 
from  aspect  sensitivity  measured  perpendicular  to  the  L-shells.  These  observa- 
tions suggest  that  the  quiet  geomagnetic  field  model  used  to  determine  aspect 
sensitivity  is  significantly  perturbed  by  the  presence  of  ionospheric  currents. 

If  the  resulting  aurora  is  an  image  of  the  bounds  of  the  auroral  electrojet  as 
suggested  by  [11,  12],  then  the  top-down  maps  provide  a picture  of  the  basic 
large-scale  structure  of  the  current  system. 

Our  theory  is  that  there  may  be  a basic  difference  in  the  magnetic  field 
distortion  induced  by  the  narrow  discrete  arcs  as  compared  to  that  induced  by 
the  large-scale  diffuse  aurora.  This  is  simply  that  the  narrow  discrete  arcs 
act  like  line  currents,  and  that  the  region  where  the  field  lines  close  around  the 
arc  is  localized  within  the  region  of  observation.  It  the  arc  is  aligned  with  the 
L-shell,  induced  field  distortion  would  occur  primarily  to  the  field  components 
which  are  transverse  to  the  direction  of  current  flow;  i.  e.  , the  L-shell  con- 
tours. This  theory  also  assumes  that  the  current  system  associated  with  the 
diffuse  aurora  resembles  a large  current  sheet. 

The  field  induced  by  this  sheet  i:  significant  in  several  ways.  First, 
the  effects  are  no  longer  localized  but  cover  a large  area,  and  the  field  lines 
do  not  encircle  the  aurora  in  the  same  fashion  as  with  the  discrete  arc.  Thus, 
within  the  reflecting  mass,  where  the  edge  effects  are  not  observed,  the  sheet 
current  associated  with  large  diffuse  auroral  forms  causes  little  perturbation 
of  the  magnetic  Z-component.  The  field  lines  from  the  sheet  will  be  generally 
aligned  perpendicular  to  the  L-shell.  The  field  strength  will  not  decrease  with 
altitude  above  and  below  the  sheet  as  in  the  case  of  the  line  current  in  which 
varies  inversely  with  range.  This  means  that  the  diffuse  aurora  creates  a 
type  of  field  distortion  which  is  fairly  uniform,  both  geographically  and  with 
altitude,  over  the  entire  viewing  region  of  the  radar.  This  type  of  field  distor- 
tion will  tend  to  shift  all  of  the  off -perpendicular  contours  of  the  radar  together 
with  no  significant  change  in  shape  or  separation  of  the  contours.  There  will 
hence  be  little  change  to  the  slope  of  the  aspect  sensitivity  curve  generated  using 
a quiescent  magnetic  field  model. 

The  more  complex  and  localized  current  systems  of  the  discrete  arcs 
can  result  in  extreme  distortion  to  the  off -perpendicular  contours  with  an 
accompanying  change  in  apparent  aspect  sensitivty  when  the  quiescent  field  model 
is  used. 

To  determine  the  feasibility  of  this  explanation,  a simplistic  model  was 
used  to  determine  the  gross  quantitative  effects  which  might  be  anticipated.  The 
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model  allows  several  line  currents  constrained  to  constant  latitude  and  altitude 
to  be  superimposed  onto  the  quiet  geomagnetic  field.  It  was  determined  that  a 
single  line  current  of  10^  to  10^  amperes  at  an  altitude  near  the  center  of  the 
arc  would  produce  variations  in  the  geomagnetic  field  similar  to  those  observed 
on  the  Fort  Churchill  magnetometer  data.  Figure  6-1  shows  a set  of  off-per- 
pendicular  contours  at  105  km  for  the  quiet  geomagnetic  field.  An  eastward 
flowing  current  of  10^  amps  at  110  km  altitude  and  53.5°  latitude  creates  a 
region  of  high  aspect  gradient  in  the  immediate  vicinity,  with  minor  perturba- 
tion north  of  the  current  and  approximately  1.5°  of  increased  aspect  angle 
immediately  south  of  the  current  merging  with  the  original  field  further  south. 
This  is  shown  in  Figure  6-2. 

If  a quiescent  field  model  were  used  to  determine  the  apparent  aspect 
sensitivity  when  the  distortion  was  similar  to  that  shown  in  Figure  6-2,  the 
slope  of  the  apparent  aspect  response  would  be  steeper  than  if  the  actual  field 
model  were  used.  This  effect  could  account  for  the  apparent  increase  in  aspect 
sensitivity  observed  in  discrete  arcs.  The  edges  of  diffuse  aurora  will  also  be 
affected  in  this  way  because,  at  these  edges,  the  current  no  longer  acts  like  a 
sheet,  and  the  field  lines  above  and  below  the  arc  connect  along  the  edge.  The 
resultant  field  distortion  acts  much  like  that  for  discrete  aurora. 


QUIESCENT  GEOMAGNETIC  ASPECT  CONTOURS  AT  105  KM  FOR  PAR 


PERTURBED  ASPECT  CONTOURS  AT  105  KM  FOR  PAR 
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A number  of  potential  problems  have  been  identified  which  could  seriously 
affect  the  determination  of  the  aspect  sensitivity  response  function.  Techniques 
have  been  described  for  elimination  or  minimization  of  many  of  the  detrimental 
effects  of  these  problems.  High  resolution  measurements  made  with  one  of  the 
most  powerful  radars  in  the  world  have  been  used  to  generate  auroral  reflectivity 
maps  and  scatter  plots  which  show  the  magnetic  aspect  response  relationship.  A 
technique  was  described  which  allowed  the  center  of  the  reflecting  arc  to  be 
accurately  estimated  so  that  only  the  most  accurate  data  would  be  selected  for 
processing. 

By  performing  pulse -by-pulse  averaging  and  by  selecting  echoes  from 
homogeneous  regions,  the  resulting  points  on  most  of  the  scatter  plots  have 
a tightly  distributed  relationship  from  which  the  "apparent"  aspect  response  can 
be  determined  with  high  confidence.  Despite  attempts  to  minimize  the  effects 
of  magnetic  field  distortion,  a wide  variation  was  found  in  the  slopes  of  the 
scatter  plots  from  different  scans  and  even  from  different  regions  in  the  same 
scan.  This  phenomena  is  apparently  related  to  the  spatial  characteristics  of 
the  sample  region,  with  a clear  cut  and  consistent  difference  between  discrete 
and  diffuse  forms.  In  the  case  of  discrete  arcs,  it  appears  that  when  the  echo 
samples  are  selected  from  a region  elongated  perpendicular  to  the  L-shells 
that  the  computed  slope  of  the  aspect  response  is  much  steeper  than  when  the 
region  is  elongated  along  the  L-shells.  This  contrast  has  not  been  observed 
for  diffuse  forms.  A theory  was  presented  which  relates  the  phenomena  to 
the  magnetic  field  induced  by  the  ionospheric  current  system.  The  results  of 
this  study  support  previous  observations  [6,  10]  that  the  apparent  aspect  sensi- 
tivity seems  to  decrease  at  greater  off -perpendicular  angles.  It  has  been  shown, 
however,  that  the  phenomena  may  result,  at  least  in  part,  from  the  superposi- 
tion of  sidelobes  and  main  beam  echo  power  causing  the  apparent  aspect  slope 
flattening  to  be  "illusory.  " 

The  findings  noted  here  clearly  indicated  that  radar  observations  of 
discrete  arcs  "appear"  much  more  aspect-sensitive  than  for  diffuse  aurora. 

It  was  also  shown  that  when  the  electrojet  current  is  confined,  as  with  the  dis- 
crete arc,  the  magnetic  field  distortion  can  increase  the  "apparent"  aspect 
sensitivity.  The  apparent  variation  in  the  magnetic  aspect  response  may  be 
merely  a manifestation  of  our  inability  to  simultaneously  measure  and  correlate 
all  of  the  spatial,  magnetic,  and  temporal  conditions  which  are  critical  to  the 
determination  of  true  aspect  response  function(s).  The  very  fact  that  aurora 
exists  means  that  the  magnetic  field  is  distorted  by  the  resulting  currents. 

If  the  relationships  suggested  here  are  valid,  and  unique  and  consistent 
aspect  response  relationships  were  developed,  then  techniques  may  be  developed 
to  use  the  differences  between  the  actual  and  apparent  magnetic  aspect  angles 
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to  determine  quantitatively  the  structure  of  the  auroral  current  systems.  This 
technique  would  be  fast  and  relatively  inexpensive  and  might  eliminate  the  need 
for  certain,  more  expensive  methods  employing  direct  measurement. 
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